Vibronic spectroscopy of an artificial molecule 
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With advanced fabrication techniques it is pos- 
sible to make nanoscale electronic structures that 
have discrete energy levels. Such structures are 
called artificial atoms because of analogy with 
true atoms. Examples of such atoms are quan- 
tum dots in semiconductor heterostructures ' 
and Josephson-junction qubits > 4 > 5 > 8 > 7 . It is also 
possible to have artificial atoms interacting with 
each other. This is an artificial molecule in the 
sense that the electronic states are analogous to 
the ones in a molecule . In this letter we present 
a different type of artificial molecule that, in ad- 
dition to electronic states, also includes the ana- 
log of nuclear vibrations in a diatomic molecule. 
Some of the earlier experiments could be inter- 
preted using this analogy, including qubits cou- 
pled to oscillators 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 and qubits 
driven by an intense field 18 ' 19 ' 20 ' 21 ' 22 . In our case 
the electronic states of the molecule are repre- 
sented by a Josephson-junction qubit, and the 
nuclear separation corresponds to the magnetic 
flux in a loop containing the qubit and an LC os- 
cillator. We probe the vibronic transitions, where 
both the electronic and vibrational states change 
simultaneously, and find that they are analogous 
to true molecules. The vibronic transitions could 
be used for sideband cooling of the oscillator, and 
we see damping up to sidebands of order 10. 

Consider two different electronic states of a diatomic 
molecule. For both states one can form a potential en- 
ergy U (r) as a function of the separation r of the two 
nuclei. In a bound state they may be approximated by 
a Morse potential, which is close to a harmonic oscil- 
lator potential near the equilibrium separation, but at 
larger distances the potential is less steep and levels off 
to allow for dissociation of the molecule. Because of cou- 
pling to the electromagnetic field, there can be transi- 
tions between the electronic levels. Associated with this, 
there is also often a change in the vibrational state of the 
molecule, and such transitions are called vibronic. The 
vibronic transitions are preferred if the two electronic 
states correspond to different equilibrium separations of 
the nuclei. The intensity of different transitions is gov- 
erned by Franck-Condon principle ' . In its classical 
form the principle says that transitions are possible be- 
tween vibrational states whose trajectories intersect in 
phase space and they are most intense between states 
having coincident turning points. 




FIG. 1: The potential energy of the artificial molecule, equa- 
tion (1). a, The total potentials (red) are formed as a sum of 
qubit energy (blue) and a harmonic vibrational potential. The 
six lowest energy eigenstates are shown for both electronic 
states. The bars represent the energy differences AE = hu 
corresponding to the drives of low and high frequency, i^lf 
and z/hf , respectively. The parameters are for the bias point 
($6,n 9 ) = (0.60<l?o, 1.16) marked by a cross in Fig. 3a. b, A 
close-up at the lowest eigenstates together with wave func- 
tions (shaded). The strongest transitions according to the 
Franck-Condon principle are shown by arrows, and are la- 
beled by the change of vibrational quanta k. 



The corresponding potentials in our artificial molecule 
are shown in Fig. 1. There are two different total poten- 
tials, shown by red lines in Fig. la, corresponding to two 
different electronic states. In general, the potentials have 
different force constants and their minima are shifted. 
The potentials are also anharmonic at higher energies. 
The only qualitative difference to real atomic potentials 
is that in our artificial molecule the leveling off of the po- 
tential, and the subsequent dissociation, is not achieved 
at the parameters we depict. Fig. lb shows a closer look 
at the potential minima, together with the lowest en- 
ergy eigenstates and the eigenfunctions. The strongest 
transitions according to the Franck-Condon principle arc 
shown by arrows. 

We have made spectroscopic studies of the artificial 
molecule. By a low frequency drive we excite the vi- 
brations and by a simultaneous high frequency drive we 
can induce the vibronic transitions. By monitoring the 
reflection of the low frequency wave we can resolve the 
vibronic levels. 

The circuit of our artificial molecule is shown in Fig. 2. 
The heart of it is a single-Cooper-pair transistor (SCPT). 
It consists of a 2 /im long Al island surrounded by two 
Josephson junctions with coupling energies Ej% w Ej2, 
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FIG. 2: SEM micrographs and circuit diagram of the artificial 
molecule, a, General view showing top electrodes of the the 
oscillator capacitors (2 large pads) with capacitances 2C each, 
b, Enlargement of the white square area in a showing the 
300 /im-long microstrip loop that forms the inductance L of 
the oscillator, c, Magnification of the white diamond in b 
showing the single Cooper pair transistor (SCPT) formed by 
two Josephson junctions with coupling energies Eji and Ej2 
and the gate electrode. Also visible is a reference SCPT (lower 
left) which enables to estimate the resistances of the junctions 
in the sample, d, The circuit diagram shows the SCPT and 
the LC oscillator connected in parallel. The resistance R = 
2.8 kfi represents the AC-dissipation of the capacitors as well 
as of other parts of the circuit. The artificial molecule is 
biased by a gate voltage V g and an external magnetic flux 
For measurement it is connected to a low frequency drive (LF) 
via a Zo = 50 transmission line and a coupling capacitor 
C c = 7 pF. The high frequency drive (HF) is used to excite 
the vibronic transitions. 



as well as capacitances C\ and G%, respectively. The 
SCPT is fabricated using standard procedure of two an- 
gle evaporation '. The SCPT is coupled in parallel to 
an LC oscillator. The inductor of the oscillator is formed 
by a superconducting on-chip microstrip loop of length 
300 /im giving L = 169 pH. The capacitance C — 196.2 
pF is formed by two parallel plate capacitors acting in se- 
ries. These are made by patterning two 1.3 mm x 1.3 mm 
Al plates on the Si3N4/Nb/Si sandwich substrate. The 
bottom plate of the capacitor is formed by a 150 nm thick 
niobium ground plane layer under the dielectric SiaN4 
layer, whose thickness is around 300 nm. 

Once fabricated, the artificial molecule has still two 
parameters that can be adjusted in the experiment. The 
dc gate voltage V g changes the gate charge en g = C g V g , 
which controls the tunneling of Cooper pairs to and from 
the island. The other parameter is the magnetic flux bias 
<!>(, through the inductor loop caused by an external mag- 
netic field. This controls the Josephson coupling Ej of 
the island as the flux shifts the relative phases of the two 
junctions and leads to complete suppresion at = $0/2 
in a symmetric SCPT. Here <& = h/2e is the flux quan- 
tum. The properties of the SCPT can be calculated in 



detail using Mathieu functions 21 ' and including the asym- 
metry of the junctions. Instead of the general analysis we 
present here a simplified treatment that is still sufficient 
to understand the main properties of the circuit. In the 
charge limit Ej < E c = e 2 /(d + C 2 + C g ), the SCPT 
forms a qubit where only two charge states are of impor- 
tance. They correspond to having one Cooper pair more 
or less on the island. In this basis the Hamiltonian of 
the SCPT is H q = \[a z E c \ - a x E JQ cos(7r<I>/<I>o)], where 
Ejo = Eji + Ej 2 , E c i = 4E C (1 - n g ) and er, are Pauli 
matrices. Diagonalizing the Hamiltonian gives two eigen- 
values E\ = ±\yjE\ + Ej cos 2 (7r<I>/<i>o), which corre- 
spond to two "electronic" states of the artificial molecule. 
The LC circuit forms a harmonic oscillator, whose Hamil- 
tonian H v = q 2 /2C+($-$ b ) 2 /2L. Taking the flux $ as 
a generalized coordinate, the charge q on the capacitor is 
the corresponding canonical momentum. Thus the total 
potential energy U of the " molecule" is 

This is plotted in Fig. 1 using parameters appropriate for 
our circuit. We see that the electronic and vibrational de- 
grees of freedom are coupled as the electronic potential 
also depends on the flux 4>. The LC oscillator potential 
vanishes at the bias point <!>{, but because of opposite 
slopes of the qubit energies, the minima of the total po- 
tentials are shifted in opposite directions > . Because 
of the different curvatures of the qubit potentials, the vi- 
brational resonance frequencies v± in the two electronic 
states are different 11 ' 13,14 . 

For measurement the molecule is connected to a trans- 
mission line. The measurement is done by sending in a 
low frequency signal typically at the frequency ^lf = 872 
MHz. This is slightly below the LC oscillator resonance 
frequency vq — 874 MHz. The reflected signal was con- 
tinuously monitored using a network analyzer set to de- 
tect magnitude and phase in a band width of 10 kHz. 
Being close to resonance, the reflected signal depends es- 
sentially on the vibrational frequency determined by the 
curvature of the qubit energy. As this varies as a function 
of the bias point (<&;,, n g ), one can use the data on the 
electronic ground state to fit the parameters of the qubit. 
We find E c /k B = 0.8 K, E J0 /k B = 3.2 K, and the asym- 
metry parameter d = (Eji — E J2 )/(E J1 + E J2 ) = 0.04. 

In order to study the vibronic transitions, a second 
drive operating at a high frequency ^hf was applied. 
Figure 3 shows the results for i/hf = 22 GHz. The re- 
flection coefficient is plotted over the (<!>;,, n g ) bias plane 
for two different low-frequency power levels Plf- The 
main observation is the yellowish or light blue fringes. 
There is one fringe starting at (<&f,,n g ) = ($0/2, 1.32). 
With increasing more fringes appear as they simul- 
taneously are curving down. One sees that the fringes 
form concentric and equally spaced ellipses around the 
point ($0/2, 1). We associate the fringe denoted by the 
dashed line with a pure electronic transition, and the 
other fringes with vibronic transitions. The concentric 
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FIG. 3: Measured spectra of the artificial molecule, a, The 
magnitude of the reflection coefficient |T| as a function of 
the bias point ($j,,n s ). The color bar gives the scale for 
r|. The dashed arch indicates the pure electronic transition 
AE/h = fuF = 22 GHz. The cross denotes the bias point of 
Fig. 1 where k — —2. The low frequency power Plf = —129 
dBm. b, the same as a but at higher Plf = —123 dBm. c, 
A cut along n g — 1.27 as indicated by the horizontal lines in 
panels a and b. d, A cut along ${, = 0.56&o- 



structure can be understood because in this region the 
energies are to a good accuracy given by Eq. (1) and, 
moreover, cos(7r<I>/<I>o) ~ — 7r(*E> / <& — §)■ This means 
that the difference in the qubit energies increases linearly 
with the distance in the vicinity of the degeneracy point 
(3>o/2, 1). In other words, the qubit bands have a conical 
shape, and the fringes correspond to approximately con- 
stant energies. The radii of the fringes are in agreement 
with the imposed frequencies and the parameters of the 
qubit, ^hf = (E+ — E^_)/h + kvo with integer k. We 
call the vibronic transitions inside of the dashed line as 
upper sidebands and outside as lower sidebands, as they 
are relative to the pure electronic transition at the same 
bias point (k > and k < 0, respectively). 

The visibility of the fringes is consistent with the 
Franck- Condon principle illustrated in Fig. 1. In par- 
ticular, at a half flux quantum (<!>;, = <l>o/2) only the 
pure electronic transition is visible. This is to be ex- 
pected since the minima of the potentials U± occur at 
the same <!> and thus only diagonal matrix elements in 
the vibrational states are nonzero. For increasing flux 
the potential minima get shifted, which allows vibronic 
transitions. The fringes appear in a crescent-shaped re- 
gion, i.e. in a stripe that is centered at the basic fringe 
(dashed line) and whose width grows with <&b (together 
with the shift of the minima of [/_ and U+). It is just 
in this region that the classical trajectories of the oscil- 
lators intersect and lead to non-vanishing transition ma- 
trix elements. The boundary of this region corresponds 



to trajectories touching each other at coincident extrema 
^cxt _ ^cxt m ^. ne £ w0 e i ec t ron i c states. Equivalently, 
for the bias points at the boundary, the oscillation just 
reaches the locus of pure electronic transitions (dashed 
line in Fig. 3). This allows to determine the oscillation 
amplitude corresponding to a given drive power Plf- 

For further analysis of the fringes we need to look more 
closely at the measurement details. In contrast to exper- 
iments with real molecules, we do not detect directly the 
absorption or emission of the vibronic quanta. Instead, 
we measure the reflection of the low-frequency wave. Two 
effects on the signal can be separated. 1) The first one is 
a usual resonance absorption taking place when the wave 
frequency is close to the vibrational frequency v±. This 
is different in the upper and lower qubit states but the 
measurements are consistent with a "motional narrow- 
ing" picture that there is only one Lorentzian line whose 
resonance frequency is an average of z/_ and v+ weighted 
by the occupations of the qubit states. This is essential 
for the visibility of the fringes near ss <& /2 as both 
the upper and lower state frequencies would be far away 
to give any noticeable absorption. 2) The second effect is 
a damping or amplification caused by the high-frequency 
drive 2 '' 28 . Consider the lower sideband fringes in Fig. 3. 
The absorption of a high-frequency quantum brings the 
system to a lower vibrational state. When the electronic 
part of the excitation relaxes by another route, this leads 
to damping of the vibrations. In a steady state this leads 
to increased absorption as the low-frequency drive has to 
supply the missing vibrational quanta. Conversely, in the 
case of upper sideband the absorption of a high-frequency 
quantum produces extra vibrational quanta leading to 
amplification of the oscillation. This could lead to a low- 
frequency reflection coefficient exceeding unity and even 
to lasing 2 ' . In our system system the losses in the vi- 
brations cause that we see this only as a decrease in the 
absorption. 

The damping/amplification effect of the high- 
frequency drive provides a natural explanation for the 
asymmetry of the upper and lower sideband fringes. In 
order to study this quantitatively, we have made numer- 
ical simulations using Bloch equations for the qubit and 
coupled to classical equations of the circuit. The simu- 
lated reflection amplitude |r| is shown in Fig. 4a. Fig. 
4b shows the energy flow from the qubit to the oscilla- 
tor. This is defined as Plf,r — -Plf, A, the low frequency 
dissipation in the oscillator resistance minus the absorp- 
tion from the LF drive in the steady state. We see that 
this is negative for lower sidebands and positive for up- 
per sidebands, in accordance with the expectation. The 
damping on the lower sidebands implies cooling, i.e. the 
microwaves lead to cooling of the oscillator, similar to 
side-band cooling in optics 2 ' . Thus we see the sideband 
cooling effect of the vibronic transitions. With a differ- 
ent measurement method, not driving the object to be 
cooled, one should get damping of thermal fluctuations. 

We comment on some additional features of the spec- 
tra in Fig. 3. There is a broad feature starting around 
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FIG. 4: Simulated properties of the artificial molecule, a, The 
reflection amplitude |T| calculated as a function of the bias 
point in the flux-charge plane. The parameters correspond 
to the measurement plotted in Fig. 3a. b, The energy flow 
from the qubit to the resonator showing damping for lower 
sideband transitions (blue) and amplification for the upper 
sideband transitions (red). The power unit is Plf = —129 
dBm. 



(0.54<j>o,l) and shifting to larger with increasing n g . 
This is the line where the vibrations in the lower elec- 
tronic state are in resonance with the low frequency ex- 
citation, v_ = i^lf- This resonance is present also in the 
absence of microwave excitation and was used to fit the 
parameters of the qubit. On the side of smaller <!>;, the 
vibronic transitions are well visible as increased absorp- 
tion because of motional narrowing since v_ < ^lf < v + ■ 
On the opposite side v_ is shifted above ^lf and the 
fringes are less clear. There are resonances also around 
(0.55<i>o, 1). They also occur in the absence of high- 



frequency drive. They are caused by Landau-Zener tun- 
neling at the point ($0/2, 1), where the electronic energy 
difference is at minimum ' . 

The data were taken at the base temperature of the 
refrigerator T — 25 mK. The high frequency power 
Phf = +15 dBm is measured at the microwave source. 
Fig. 3 shows the absolute value of the reflection coefficient 
defined by V^ ut = TV in . In addition to several stages of 
band pass filtering by Minicircuits VLX-series filters, we 
used two circulators at 25 mK to prevent the back-action 
noise of our cryogenic low-noise amplifier from reaching 
the qubit. The attenuation of the rf- signal lines were 
measured at room temperature. Cooling down to 25 mK 
reduced the attenuation by 6 dB. We checked that the 
AC Stark shift at large drive powers agreed with the ex- 
citation amplitudes on the measurement drive. 

The simulation uses the circuit diagram of Fig. 2 with 
Ti = 40 ns, T2 = 2 ns and other parameters mentioned 
in the text. In producing Fig. 4 the cone approximation 
mentioned in the text is used as the difference to full 
Mathieu bands is negligible. 
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